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We discuss possible experimental signatures of forming a Quark-Gluon plasma in 
high energy nuclear collisions. In first order phase transitions such as the chiral 
^ ' symmetry restoration supercooling may lead to density fluctuations such as droplet 

^■f^ ' formation and hot spots. These will lead to rapidity fluctuations event-by-event 

and triggering on such fluctuations will have distinct effects for HBT radii. 

o 

a^ 

I 1 Introduction 

OO 

, Results from the RHIC collider next year are eagerly awaited. The hope is 

to observe the phase transition to quark-gluon plasma, the chirally restored 
hadronic matter and/or deconfinement. This may be by distinct signals of J/^ 

, , suppression, strangeness enhancement, rj' enhancement, enhanced temperature 

O , and multiplicity fluctuationsu, etc. Other signals are plateau's in temperatures, 

transverse flow or other collective quantities as function of centrality, transverse 
energy or multiplicity. We will here suggest other signals of a first order phase 
transition namely rapidity fluctuations and associated variations in HBT radii. 



5— ( ' 2 Supercooling and Droplet formation 



Lattice calculations suggest that QCD has a first order transitioirat zero chemi- 
cal potential provided that the stranjge quark is sufficiently lightH. Recent work 
using random matrix theory (RMT)El suggests an effective thermodynamic po- 
tential of the form 

n{cf>)/Nf = 02 - iln{[</.2 _ ^ ^y)2j . [^2 _ _ ^y)2j| 

The value oi (j> ^ (V'V') at the minima of is related to the expectation 

value of the quark condensate which is the order parameter for chiral symmetry 
breaking. Minimization of Eq. (Q) leads to a fifth order polynomial equation for 
(j) which is identical in form to the results of Landau-Ginzberg theory using a 4>^ 
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potential. One solution to this equation corresponds to the restored symmetry 
phase with = 0. This model predicts a second order transition for fi — 
at a temperature, Tc, which is generally agreed to be approximately 140 Mev. 
For T = 0, a first order transition occurs at some ^o- Since the phases in 
which chiral symmetry is broken and restored must be separated in the (/i, T) 
plane by an unbroken line of phase transitions, this implies the existence of a 
tricritical point. In this model, the tricritical point is at 



Further, the discontinuity in the density at (T = 0,/io) can be estimated 
as An w 2.5no where no is the equilibrium density of nuclear matter. An 
inevitable consequence of such a phase diagram is the existence of spinodal 
lines within which ^l{(f)) has three local minima. In this region, nuclear matter 
can be either superheated or super cooled. 

The qualitative features of this phase diagram lead to an interesting sce- 
nario for relativistic heavy ion collisions in which matter is compressed and 
heated. Some or all of this matter will undergo chiral restoration. If the 
subsequent expansion is sufficiently rapid, matter will pass the phase coexis- 
tence curve with little effect since the restored symmetry phase with (f) — 
will remain a local minimum of Q{(t>). This suggests the possible formation of 
"droplets" of supercooled chiral symmetric matter with relatively high baryon 
and energy densities in a background of low density broken symmetry mat- 
ter. These droplets can persist until the system reachs the spinodal line and 
then return rapidly to the now- unique broken symmetry minimum of i~t{(f)). 
The schematic nature of Eq. (|^) makes it impossible to make quantitative pre- 
dictions about the properties of such droplets. (The effective thermodynamic 
potential can well have additional terms which are independent of ((>.) However, 
a large mismatch in density and energy density seems to be a robust predic- 
tion. Note that these "droplets" are kineniatically distinct from the "hot spots" 
usually associated with jets and minijetsu. Jets are characterized by a strong 
directional orientation and a high transverse momentum not expected for the 
present droplets. 

Such density fluctuations can occur only for first order transitions. In 
the RMT model, these will only be the case for sufficiently high net baryon 
densities. In ultrarelativistic nuclear collisions at mid rapidity, the baryon 
density is low as in the early universe and, according to RMT, the transition is 
of second order. One should then investigate nuclear fragmentation regions or 
go to lower collision energies (e.g., AGS energies) to observe droplet formation. 
However, most lattice simulations do find a first order transition at zero baryon 
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density. 



3 Rapidity Fluctuations 

If the transitions is first order, matter may supercool and subsequently create 
fluctuations in a number of quantities. Density fluctuations in the form of 
hot spots or droplets of dense matter with hadronic gas in between is a likely 
outcome. We shall refer to these regions of dense and hot matter in space-time 
as well as in momentum space as droplets. If we assume that hadrons emerge 
as a Boltzmann distribution with temperature T from each droplet and ignore 
transverse flow, the resulting particle distribution is 



dy(Pp±_ 



Here, y is the particle rapidity and p± its transverse momentum, fi is the 
number of particles hadronizing from each droplet i, and 

2 ti- Zi 2 l-Vi 

is the rapidity of droplet i. 

When m±/T 3> 1, we can approximate cosh(j/ — rji) ~ 1 + ^{y — rji)^ in 
Eq.(H). The Boltzmann factor determines the width of the droplet rapidity 
distribution as ^ yjT /m±. The rapidity distribution will display fluctuations 
in rapidity event by event when the droplets are separated by rapidities larger 
than \T]i — rij\^y^T/m±. If they are evenly distributed by smaller rapidity 
differences, the resulting rapidity distribution (|^) will appear flat. 

The droplets are separated in rapidity by \r]i — r]j\ ^ Az/tq, where Az is 
the correlation length in the dense and hot mixed phase and tq is the invariant 
time after collision at which the droplets form. Assuming that Az ~ Ifm — a 
typical hadronic scale — and that the droplets form very early ro^lfm/c, we 
find that indeed \rii — rij\}ty^T/m± even for the light pions. 

If strong transverse flow is present in the source, the droplets may also 
move in a transverse direction. In that case the distribution in pj_ may be 
non-thermal and azimuthally asymmetric. 
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Figure 1; Droplet formation (top) and corresponding rapidity distribution (bottom). 

4 HBT and Correlation Functions 

In this section we briefly review the correlation function analysis of Bose- 
Einstein interference. For details we refer, to references on aiellar intensity 
interferometryEl, pion interferometry in ppcl and AA collisions LtlIij. 

For a source of size R we consider two particles emitted a distance ~ R 
apart with relative momentum q = (ki — and average momentum, K = 
(ki +k2)/2. Typical heavy ion sources in nuclear collisions are of size i? ~ 5 fm, 
so that interference occurs predominantly when q^fi/R ^ 40 MeV/c. Since 
typical particle momenta are ki/tK ^ 300 MeV/c, the interfering particles 
travel almost parallel, i.e., ki k2 — K ^ q. The correlation function due 
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to Bose-Eins-tein interference of identical particles from an incoherent source 
is (see, e.g.,B) 

/ d^x S{x, K) e*9^ ^ 



C2(q,K) = l ± 



(4) 



Jd^x S{x,K) 

where S{x, K) is the source distribution function describing the phase space 
density of the emitting source. The +/— refers to boson/fermions respectively. 

Experimentally the correlation functions for identical mesons (tt^tt^, 
K^K^, etc.) are often parametrized by the gaussian form 

C2iqs,qo, m) = 1 + Aexp [-qIrI - qlRl - qfRf ~ iq^qiRli] . (5) 

Here, q = ki — k2 — {qs,qo,qi) is the relative momentum between the two par- 
ticles and Ri, i = s,o,l the corresponding sideward, outward and longitudinal 
HBT radii respectively. We have suppressed the K dependence. We will em- 
ploy the standard geometry, where the longitudinal direction is along the beam 
axis, the outward direction is along K, and the sideward axis is perpendicular 
to these. Usually, each pair of mesons is lorentz boosted longitudinal to the 
system in which their rapidity vanishes, Y = 0. Their average momentum K is 
then perpendicular to the beam axis and is chosen as the outward direction. In 
this system the pair velocity (3k='^/ Ek points in the outward direction with 
/3o = p±/rn± where m± = ^ + p'^ is the transverse mass. Often models 
assume a boost invariant source and boost to the cms where y = where a 
number of boost factors disappear and the out-longitudinal coupling Roi van- 
ishes. The reduction factor A in Eq.(^) may be duCj-to long lived resonances 
eH13, coherence effects, incorrect Gamow corrections or other effects. It is 
found to be A ~ 0.5 for pions and A ^ 0.9 for kaons. 



5 HBT for Droplets 

Droplets lead to spatial fluctuations in density which can be probed by correla- 
tions between identical particlesu. For simplicity we parametrize these droplets 
by spatial and temporal gaussians of size Rd and duration of emission td- The 
distribution of particles in space and time — usually referred to as the source 
— for droplets situated at Xi = (r^, i, ) is 

S{x,K) - ^^(a;,,is:)exp 

i 

where S{xi^ K) is the distribution of sources. Normalizations of S cancel when 
calculating correlation functions, Eq.(^). Other effects as transverse flow, opac- 
ities llHl, resonances t3, etc. can be included but we shall ignore them here for 



(r - r,y {t - t,Y 



2*2 



(6) 
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simplicity. Likewise, the extension to droplets of different size and duration of 
emission is straight forward but unnecessary for our purpose. We expect the 
scale of S is the nuclear overlap size Ra of order several ferrnVs whereas the 
droplet size Rd is smaller — of order a few /ermi's. 
From and (H) we obtain 

C2{q) ^ 1 + cxphq^i?^ - (q . pftl] \~S{q, K)\'' , (7) 

where 

S{q, K) = Y1 S{x^,K)e^''-^' I J2 Six., K) (8) 

i i 

is the Fourier transform of the distribution of sources and (3 — 'K./Kq is the 
average velocity of the pair. 

If there is only a single droplet S{q, K) = 1 and the correlation function is 
simply given by the droplet gaussian. If there are many droplets, the sum can 
be replaced by an integral. 

5.1 Non-expanding Sources 

For non-expanding sources we can assume that the particle momentum distri- 
bution is the same for all droplets, i.e., independent of K. 

For many droplets we assume that the droplets are distributed by a gaus- 
sian 

~S{x,, K) ^ exp(-r2/2i?i - ty2R\) , (9) 

where Ia is the spread in droplet formation times and Ra is the size of the 
nuclear overlap zone. It increases with decreasing impact parameter - roughly 
as RA{b) ~ .JRa{QY - 62, where i?A(0) ~ l.2A^/'^fm is the nuclear size. 
The resulting correlation function becomes 

C2{q) ^ 1 + exp[-q2(i?i + _ (q . p)\t\ + t^)] . (10) 

Since Ra and cIa are typically of nuclear scales ~ 5 — 10 fm whereas we expect 
Rd ''^ 1 fm, the droplets are simply "drowned" in the background of the other 
droplets. 

In the case of only two droplets (of same size) the resulting correlation 
function is 

C2iq) ^ 1 + exp[-q2i?2 _ (q . pff}] cos(iq • (xi - X2)) . (11) 

However, as the Xi's differ from event to event the oscillation in (^l]) will, 
when summing over events, result in a Fourier transform over the distribution 
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of Xi^s in different event and the end result will be similar to (p^. The number 
of pairs in a single event at RHIC energies is probably not large enough to 
do event-by-event HBT. The oscillation will also be smeared by resonances, 
Coulomb effects, a hadronic background, etc. 

We conclude that it will be very difhcult to see the individual droplets 
through HBT when the sources do not expand. 



5.2 Expanding Sources 

At RHIC energies the collision regions are rapidly expanding particularly in the 
longitudinal direction. Consequently the droplets may have different expansion 
velocities and rapidities and the distribution S{xi,K) will depend on K and 
this difference will now be exploited. 

When the droplets are separated in rapidity such that \m — r]j\^y/T/m±, 
only particles within the same droplet contribute to the correlation function. If 
we transform to the droplet center-of-mass system, i.e. rji — 0, the correlation 



function is then given by ( 10) with S — 1. In terms of qs, qo, qi it can be written 

C2{q) - 1 + exp [-qlRl - qliRj + ^14) - q^R^ + p^tl) - 2qoqiMltl] , 

(12) 

where Po = p±/rn±cosh{Y) and (3i = tanh(y), are the transverse (outward) 
and longitudinal velocities respectively of the pair. 

When the droplets overlap in rapidity, we need to consider the distribution 
of droplets in more detail. At ultrarelativistic energies the rapidity distribution 
is expected to display approximate Bjorken scaling, i.e., the rapidity distribu- 
tion is approximately given by (§) where the droplet rapidities rji are more or 
less evenly distributed between target and projectile rapidities. Parametrizing 
the transverse and temporal distribution as gaussians, we arrive at the droplet 
distribution 



S{xi, K) ^ exp 



cosh(y — r7i) ^ — ^ 77"!-^ 

T ^ ' ' 2R\ 2t\ 



(13) 



where — \/t'f — zf is the invariant time and Tf the average freeze-out time. 

The resulting correlation function becomes (in the system Y — 0, where 
Pi = tanhF = and f3o = p±/m±) 

C2{q) = 1 + expi-q^RA + Rl) - qliRl + Rl + PliA + tl)) 

-lfir'f:^+R',)]. (14) 
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R,(b) [fm] 

Figure 2: The HBT radii as function of nuclear overlap Ra which is proportional to centrality, 
dN/dy or E^. Droplet formation is assumed to set from semicentral collisions. A long mixed 
phase will then lead to large outward HBT radius Ro whereas triggering on large rapidity 
fluctuations corresponding to small droplets leads to smaller longitudinal Ri and sideward 
Rs HBT radii. 

We observe that the larger nuclear size Ra dominates the smaller droplet size 
Rd when droplets overlap. 

6 Centrality dependence of HBT radii 

We can now study the consequences of forming droplets at RHIC energies. The 
onset of large rapidity fluctuations is one signal but it should be accompagnied 
by the following behavior of the HBT radii. In peripheral collisions, where 
the nuclear overlap and stopping is small, we do not expect suSicient energy 
densities to build op and form droplets. Thus the HBT radii should simply 
grow with the size of the geometrical nuclear overlap, centrality, dN/dy and 
transverse energy Et- At SPS energies the HBT radii are indeed found to scale 
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approximately with the geometrical sizes of the colliding systems as given by 

If energy densities achieved in RHIC collisions are sufficient to form droplets 
in central collisions, the HBT radii will deviate from the geometrical overlap if 
triggered on large rapidity fluctuations (see Fig. 2) as given by Eg. (|l2[) . Com- 
paring the theoretical predictions of Eqs. (^2|) and (|lj) with the experimentally 
measured HBT radii of Eq. (||) we see that the sideward and longitudinal HBT 
radii decrease from the nuclear size tA to the droplet sizes Rd, td- Thus at 
a certain semi-centrality, where energy densities achieved in nuclear collisions 
start becoming large enough to create droplets, the sideward and longitudinal 
HBT radii should bend over and start decreasing with centrality. The outward 
HBT radius may behave differently depending on the duration of emission. 
The droplets may emit hadrons for a long time, as is the case for a long lived 
mixed phase-Xreferred to as the "burning log" ) . In the hydrodynamic calcula- 
tion of Ref. the duration of emission and consequently the outward HBT 
radius increase drastically up to five times larger than the transverse size of 
the system, i.e. Ra- This scenario is indicated in Fig. 2. 

Besides droplets we may expect some hadronic background. It is straight 
forward to include such one in the correlation function. As its spatial extend 
is expected to be on the scale ~ Ra, it will reduce the correlation function at 
small q ~ Ti/Ra- However, the droplets will still lead to correlations at large 
relative momenta of order q ^ h/Rd- The large q correlations are suppressed by 
the square of the fraction of pions emerging from droplets at a given rapidity. 

7 Summary 

If first order transitions occur in high energy nuclear collisions, density fluctua- 
tions are expected which may show up in rapidity fluctuations event-by-event. 
HBT interferometry adds an important space-time picture to the purely mo- 
mentum space information one gets from single particle spectra. Triggering on 
rapidity fluctuations, the HBT radii may display a curious behavior. The out- 
ward HBT radius Rg may increase drastically with centrality due to a long lived 
mixed phase. The longitudinal Ri and sideward Rs HBT radii will, however, 
saturate and decrease for the very central collisions because a large rapidity 
fluctuation signals a hot and dense droplet of small size. 

The predicted behavior for the sideward and longitudinal HBT radii is 
opposite to that predicted in cascade and hydrodynamic calculations. It would 
be a clean signal of a first order phase transition in nuclear collisions. 
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